High Order Coherent Control Sequences of Finite- Width Pulses 
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The performance of sequences of designed pulses of finite length r is analyzed for a bath of spins 
and it is compared with that of sequences of ideal, instantaneous pulses. The degree of the design 
of the pulse strongly affects the performance of the sequences. Non-equidistant, adapted sequences 
of pulses, which equal instantaneous ones up to C(r 3 ), outperform equidistant or concatenated 
sequences. Moreover, they do so at low energy cost which grows only logarithmically with the 
number of pulses, in contrast to standard pulses with linear growth. 
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The rapid evolution of the field of quantum science and 
quantum information demands robust quantum control 
techniques in the presence of environmental noise. To dy- 
namically generate systems essentially free from decoher- 
ence has now become a focus of the research of quantum 
control. This suppression of decoherence is an important 
requisite in quantum information processing [l[, for ex- 
ample for the realization of a quantum computer, in nu- 
clear magnetic resonance (NMR) , for high accuracy mea- 
surements [2] or in magnetic resonance imaging (MRI) 
to mention only a few. 

In this work we focus on quantum control by short 
pulses of finite length. It is beyond our scope to discuss 
continuous quantum control, see for instance Ref. H it 
is on a discovery in NMR, the Hahn spin echo [H, that 
the pulsed-control methods are based. The original tech- 
nique makes use of an electromagnetic pulse in order to 
rotate the spin and to refocus it along a desired direction. 
Dynamical decoupling (DD) 0,0] iterates the single pulse 
in a sequence of pulses such that the coupling between 
the spin and its environment is averaged to zero. Among 
the "open-loop" pulse-control techniques, the dynamical 
decoupling is one of the most promising protocols for pro- 
longing the coherence time of a spin (qubit) coupled to 
an environment. No detailed, quantitative knowledge of 
the decohering environment is required. 

The sequences come in a large variety. We distinguish 
equidistant and non-equidistant sequences. In the first 
category we recall the iterated Carr-Purccll-Mciboom- 
Gill (CPMG) sequence [1,0, where the pulses are reg- 
ularly separated (apart from the very first and the very 
last one). To the second category belong for instance 
the universal Uhrig DD (UDD) sequence [loWl^. the 
Locally Optimized Dynamical Decoupling (LODD) fl3| . 
the Optimized Noise Filtration by Dynamical Decoupling 
(OFDD) and the Bandwidth-Adapted Dynamical 
Decoupling (BADD) [l5[ for pure dephasing models and 
the concatenated DD (CDD) [llor UDD (CUDD) [l7| 
or the quadratic UDD (QDD) [18( for models with de- 
phasing and relaxation. 

The design of the DD schemes relies originally on the 



assumption that the pulses are arbitrarily strong and in- 
stantaneous though the effects of pulses of finite length 
were known to matter 1^-23|. But the pulses used in 
laboratories always have a bounded, finite amplitude so 
that they have a finite duration r. Even if sequences 
like CPMG and UDD have already been implemented in 
experiments with very good results 0, [l3 , Q£ 
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that pulses have a finite duration appears often as a nui- 
sance deteriorating the s upp ression of decoherence, see 
for instance Refs. |2fJ andl22l 

It is of great practical relevance to which extent the 
length of a pulse affects the performance of a sequence 
such as UDD or CPMG of given duration T. How should 
one choose the location, the duration (or the amplitude), 
and the shape 2l|, [23|, |25[ of the bounded pulse in order to 
minimize the errors due to its finite duration if it replaces 
the ideal, instantaneous pulses in a certain sequence? 

Here we report for the first time numerical evidence 
of how sequences of realistic pulses of finite width must 
be designed in order to achieve the same perturbative 
suppression of dephasing as the corresponding ideal se- 
quence. We compare various known sequences 2(il. 26l 27 1 
and numerically analyze their performance for a spin cou- 
pled to a bath of spins. To obtain an experimentally rel- 
evant comparison all pulses are designed in such a way 
that the largest amplitude appearing in each sequence is 
the same [28j . 

The Model. We consider the pure dephasing Hamil- 
tonian H = l g <E> Bq + o z ® B z that determines the free 
evolution of the system between two consecutive pulses 
by Uf ICC (t) = cxp{—itH}. The operators Bq and B z act 
on the bath only, while the identity l q and the Pauli ma- 
trix a z act on the qubit represented by a spin 1/2. For 



simplicity we identify henceforth l q $ 
bath consists of M spins with i € {1, . 



H = u h B + ctW J2^ a 



Bo and Bo. 
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No drift term oc a z of the qubit is included because we 
work in the rotating reference frame. Explicitly we an- 
alyze two cases, see also Fig. [TJ (i) A spin chain with 



Bo = E 



M 



CT W . ^+i) a, = A and N s = 1. (ii) A 
central spin model [ll|, [22-[i3 characterized by a dipolar 
coupling B = Y,f<i=i (aaPvP-ffV-ff®) with 
Aj (M - 1) = A(2i - M - 1) and N s = M. The rapidity 
of the dynamics of the bath is given by Wb •= ol\ with a 
a dimensionlcss constant. 
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FIG. 1: Left: Spin bath models under study with a qubit 
(square) coupled (dashed lines) to M bath spins (dots) inter- 
acting among themselves (solid lines). Right: Upper panel: 
1 st order SCORPSE tt pulse [H; middle panel: 2 nd order tt 
pulse; lower panel: 2 nd order 2tt pulse. Both 2 nd order pulses 
are found by solving the conditions derived in Ref. [25|. Am- 
plitudes and switching instants are available upon request. 




FIG. 2: (Color online) Upper panel: Sequences of type (i) 
(tj-CPMG, TjCDD or Tj UDD) are sketched. Lower panel: Se- 
quences of type (ii) are shown (tj RUDD). Only the maximum 
amplitude and the pulse duration are shown, but no details 
of the pulse shapes. The 7r-pulses are depicted by filled (red) 
blocks, the initial and final 27r-pulses by open (blue) blocks. 
The first tt pulse of TjRUDD and all pulses of type (i) se- 
quences have the same amplitude a max to ensure experimen- 
tally relevant comparability. The instants t% are given in Eqs. 
[5] and [31 the start and end points tf in Eq. Q. 



pure dcphasing. The CDD sequence of level k is defined 
by the recursion 

CDD fc+1 (T) = CDD fe (T/2) o h, o CDD fc (T/2) , (3a) 



The control Hamiltonian is given by H c (t) = cr®v(t). 
We consider piecewise constant pulses shown in Fig. [T] 
During each pulse of total length t^' the qubit evolves 
under the simultaneous action of the system and of 

the control Hamiltonian U p = Texp{— i J^° +T (H + 
H c (t))dt} where T stands for standard time ordering. 
The evolution operator of the total sequence from t = 
to t = T is denoted by R. 

The Sequences. Two types of sequences are studied, 
see also Fig. [2J (i) The durations r^' = t* of the pulses 
is constant throughout the sequence and it is kept con- 
stant on variation of T. These sequences are denoted by 
Tj CPMG, TjCDD and TjUDD, because they reproduce 
the ideal CPMG, CDD and UDD sequences for r -> 0. 
The subscript j stands for properties of the pulses as ex- 
plained below, (ii) The durations are varied along 
the sequence, i.e., they depend on i. But they shall not 
depend on T other than that the sum of all pulse du- 
rations cannot exceed T, i.e., T > T p := Ei T • The 
corresponding sequences are denoted by Tj RUDD. 

The sequences of type (i) are made of N tt pulses whose 
width is t* . The center of the i-th pulse is given by 



UDD 



CPMG 



T (2i — l)/(2iV), 



(2a) 



T sin 2 (tt i/(2(./V + 1))) , (2b) 

for the TjCPMG and r^UDD sequence, respectively. 
We use the simplified version of CDD designed only for 



CDD fc+1 (T) = CDD fe (T/2) o CDD fc (T/2) , (3b) 

where (13"a|) holds for k even and (|3b[) for k odd; o stands 
for concatenation and T\ v for the operator of a pulse of 
angle ip. The zero-level CDDo(T) is free evolution with- 
out pulses. 

The subscript j in Tj refers to the order of the 
pulses, i.e., its time evolution operator fulfills U p = 
exp{— irSojn^ + 0(r J+1 ). We restrict our study here 
to explicit pulses with j = 0, 1, 2, see Fig. [TJ which fulfill 
the conditions derived in Ref. |25| . A recursion for gen- 
eral j is given in Ref. [3ol The th order pulse is simply 
rectangular; the other pulses used are depicted in Fig. [T] 

The sequences of type (ii) are similar to the TjUDD 
sequences in that they are based on pulses of order j. 
The crucial difference is that pulses are not constant in 
leng th. They are defined according to our previous work 
[271 ] by a start instant t~ and a stop instant given by 



t± := Tsin 2 



2{N+ 1) 



± 



%{T) 



(4) 



The above relation results naturally from the requirement 
that the effective switching function of the sequence ex- 
pressed in 9 £ [0,7r] according to t = Tsin 2 ((9/2) is an- 
tiperiodic [22j- This antiperiodicity ensures that the to- 
tal sequence suppresses the decohering terms oc <r z in the 
time evolution [12[. The duration of the pulses in time 
t^ = tf — fc~ yielding 



= Tsin(7r i/(N+ l))sin(0 p ) 



(5) 



is determined by the parameter O p (T). It acquires a de- 
pendence on T if we require r* := r' 1 ' to be constant 
upon varying T. Note that 9 p = n/(2(N + 1)) refers to 
back-to-back pulses without any free evolution between 
them, see below. 

The antiperiodicity of the switching function is the ba- 
sis for the suppression of dephasing in high order [T2I l36j . 
In order to guarantee this antiperiodicity, it is required 
to insert an initial and a final pulse which represent the 
identity U p = cxp{— irBo} + Ofr^ 1 ). For instance, it 
may be a zero tt or a 2tt pulse [27| . The initial pulse starts 
at t^ = and stops at t^ = Tsin 2 (0 p /2) while the fi- 
nal one starts at t~^ +1 = Tsin 2 [(n — p )/2] and stops at 
fjy-ij = T. These pulses are indicated by open boxes in 

Fig. m 

In the sequel, we compare the various sequences always 
with the same t* because the shortest accessible pulse 
duration of a 7r pulse, corresponding to the maximum 
amp litude, represents a crucial experimental constraint 
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. Only the very short boundary 2ir pulses in the 
TjRUDD are treated separately. But their importance is 
assessed by considering TjRUDD with and without the 
boundary 2ir pulses. We stress that due to the variable 
duration of the pulses according to (|4I5|) in the RUDD 
sequence most of the pulses are much longer than r* . 

The Partial Frobenius (A p p) Distance defines the 
distance between the ideal evolution of the initial state 
of the qubit due to the pulses and its evolution includ- 
ing the interaction with the bath and the application of 
the sequence 37[. For each axis of rotation 7 = {x, y, z} 
we define a difference of density matrices of the qubit by 



(7) 



tr B 



(7) 



Jn) 



where p^ 



Rp^ft. 



The 

partial trace over the bath is denoted by tre- Given a 
factorized initial state p^ := I7) (7I ® 1b the density ma- 
trix pij '■— (?x Po'^x i s trie ideally evolved po subject 
only to ideal pulses without any bath interaction. The 
distance A p p measures the difference between the real 
evolution and the ideal one reading 



A 2 



1 \ - 



(0) 



■y=x,y,z 



Numerical Simulation. We compute the performance 
of sequences of pulses of finite duration for the systems 
in (flj shown in Fig. [TJ We choose the minimum duration 
r* < minijT^} and a minimum value of T such that 
T > J^i T ^ ' see captions for values. Sequences with 
N = 10 pulses are considered because this number allows 
us to consider the CDD sequence as well; it corresponds 
to the concatenation level k = 4, cf. Eq. ©. The results 
are shown in Figs. [3] and Hta) for the spin chain model 
and in Figs.[4jb) and Fig. [5] for the central spin model. 

First, we consider the influence of the topology and 
the size of the spin bath. In Fig. [3f c) data for the spin 
chain is shown for M = 3 (open symbols) and data for 
AI = 8 (filled symbols) fits in perfectly. This indicates 
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FIG. 3: (Color online) Distance A p p vs. the duration T of 
sequences of pulses with zero or finite width for a = 10 for 
spin chains. All open symbols refer to M — 3 bath spins; the 
filled symbols in panel (c) to M = 8 bath spins. The finite- 
width pulses have minimum width r* = 1.086 • 10~ 3 /A and 
T > 0.09/A; Panel (a): rectangular th order pulses; Panel 
(b): SCORPSE 1 st order pulses [H]; Panel ( c ) : 2 nd order 
pulses shown in Fig. [T] To highlight power-law behavior the 
dashed lines are included: The UDD curve scales as t 2 ( n+1 \ 



that the size effect is very small in the regime of interest. 
The topology of the spin bath has a certain impact, but 
only on the quantitative level, not on the qualitative one 
as can be seen comparing Fig. [3](c) with Fig. SJb). The 
results for the central spin model with M = 8 bath spins 
are qualitatively identical to the ones for the spin chain 
except for a heuristic factor k ~ 2.3 in T. The latter 
can easily be understood in the sense of an effectively 
stronger coupling between qubit and bath for the central 
spin model than for the spin chain for the same value A 
because there are more couplings Ai oc A between qubit 
and bath spins. 

Second, we study the influence of the sequences on 
the performance. Thus we consider long sequence dura- 
tions T. In this regime the pulse errors are unimportant 
and pulse shaping plays only a minor role. This fact is 
perfectly understandable because for given r* the limit 
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FIG. 4: (Color online) Same as in Fig. [3jc). Panel (a): spin 
chain with M = 3, a = 100 and r* = 1.086 ■ 10~ 4 /A, T > 
0.009/A. Panel (b): central spin model with M = 8, a = 10 
and t* = 0.0004828/A, T > 0.04/A. The result is identical to 
Fig. [31c) except for a shift by the factor k « 2.3 in T. 



T — > oo implies that r* /T vanishes. In the formalism 
of filter functions [ltj EH, [38l - l40j this can easily be seen. 
The signal s(T) = exp(— 2x(T)) is determined by the 
frequency integral 



X(T) 



1 , 1* 



(7) 



where F{uT) is the hltcr function. For pulses of duration 
t^i centered at instants SjT it is given by 



N 



F(z)= ! + (-!) 



N+l e -iz 



cos I 



2T 



, (8) 



where we use z := wT for brevity. This equation is valid 
if the coupling between qubit and bath is effectively zero 
during the pulse. For artifical noise this can be realized 
experimentally fl3j while for generic systems the pulse 
design has to approximate this situation [25, 27 1. Clearly, 
for larger and larger T the influence of the finite pulse 
durations t^> decreases more and more. 

The scaling of A p p with T for UDD with ideal pulses 

is also remarkable. For UDD, d ig \b^ | + \b--\ where 

6-1 and b depend on B$, B z and on the initial density 

matrix p [36(. In particular, cc T 2 ( N + 1 ) an d its pref- 

actor is even in B z while 6_| oc j 1 ( Ar + 1 ) with a prefactor 

odd in B z . If both are present one has the generic result 
d = 0(T N+1 ). But if the Hamiltonian is symmetric un- 
der global spin flip a z <H> — <r 2 , realized, e.g., by a 7r rota- 
tion about total a Xl it follows that b^ = — b^ =0 due 



to its oddness in B z such that we obtain d = 0(T 2(JV+1) ) 
which is better than gcnerically expected. Hence on the 
one hand, the generic behavior of dynamic decoupling 
can only be seen for systems without symmetry. On the 
other hand, we stick here to the Hamiltonian ([T]) because 
it is of the kind occurring mostly in experiment[2 , 29], 41 1 . 

Fig. (He) with a = 10 and Fig.HJa) with a = 100 differ 
in the rapidity of the bath dynamics which is faster for 
larger a. Clearly, the decoherence sets in earlier if the 
bath is faster because the switching by the pulses is rela- 
tively slower. This is no contradiction to the basic idea of 
motional narrowing stating that a very fast bath implies 
longer coherence times because the fast bath dynamics 
reduces its influence on the qubit due to averaging. But 
previous results, e.g., Fig. 3 in Ref. Il8l. show that for this 
effect to take place a should exceed 10 6 . 

We do not consider data for smaller a J$ 1 here because 
it is the our present scope to show how the detrimental 
effect of finite pulse duration can be compensated. But 
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a previous study on single pulses, see Fig. 7 in Ref. 
revealed that effects of the finite duration of the pulses 
become noticeable only for a > 1. 

Third, we consider the large regime of shorter dura- 
tions T where A p p is dominated by the properties of 
the pulses. Naturally, this effect is most prominent for 
the uncorrected rectangular pulses of order. In Fig. 
(3ja) the distance d is significantly larger for pulses of 
finite width (symbols) than for the ideal ones (lines). 
The RUDD sequence performs worse than the other se- 
quences. This is not surprising since it is based on the 
assumption that the pulse is designed such that there is 
none or no significant coupling between qubit and bath 
during the pulse. A rectangular pulse realizes this as- 
sumption only in order r* . 

Hence it is clear that the level for A p p which can be 
reached for small values of T is lower for the 1 st or- 
der pulses (panel (b)) and even lower for the 2 nd order 
pulses (panel (c)). This fact illustrates nicely that the 
optimization of pulses is indeed an important ingredi- 
ent in enhancing the performance of dynamic decoupling 
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The key observation is that the Tj-RUDD becomes the 
best performing for j = 2. For j = and j = 1 the 
TjUDD sequence turned out to be more advantageous. 
We conclude that the pulses need to be sufficiently well 
designed in order that the underlying idea of the RUDD 
sequence [13] really pays. In Fig. OJc) the gain using 
RUDD instead of UDD is about two orders of magnitude. 
Such improvements are to be expected in the regime 
where the performance of the sequences is dominated by 
the pulse errors. 

We emphasize that the fact that RUDD performs bet- 
ter than UDD or any other generic sequence of pulses 
of constant duration is quite remarkable because most of 
the pulses in the RUDD sequence are much longer than 
t* . The sum T p of the lengths of all N pulses is T p = Nt* 
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for a generic sequence while it is 

T p = T*cot(7r/(2(iV + l))/sin(7r/(JV + 1)) (9a) 
w t*2{N + 1) 2 /tt 2 for N large (9b) 



for the RUDD sequence according to Eqs. (|4l5p . One may 
prefer to consider the total energy necessary to realize 
the sequence Q. The energy required for a given pulse 
is proportional to 1/r. Hence the total energy E p is 
given for the UDD sequence by E p = AN/t* where A 
is a constant depending on the shape of the pulse. Note 
the linear divergence in N. In contrast, for the RUDD 
sequence one obtains 



Asin(7r/(iV + 1)) 



N 



1 



-sin(^/(7V + l)) 



(10a) 



« (2A/r*) In [2 (AT + 1)/tt] for AT large (10b) 

which diverges only logarithmically in AT. Thus, given 
a minimum pulse duration r* it is much less costly in 
energy to reach long coherence times by applying RUDD 
than by any generic sequence with pulses of constant r* . 

In view of the above observations, it remains to clarify 
why the RUDD works better than the other sequences, 
but only for higher order pulses. According to the ana- 
lytic foundation of RUDD [22| , its advantage over other 
sequences with shaped pulses consists in the vanishing of 
mixed terms in T and r*. For instance, an ideal UDDjy 
scales generically like T N+1 and the TjUDDjv of N finite- 
width pulses certainly has errors scaling like T N+1 and 
(t*) j+1 . But one cannot exclude the occurrence of terms 
such as TV*, TV, or T(t*) 2 . They result from the in- 
terplay between the finite duration of the pulses and the 
sequence. It is crucial that this is different for TjRUDDjv. 
There the finite duration is fully taken into account in 
the design of the sequence 27| . Hence the errors of the 
Tj RUDD N are of the order T N+1 and (r*) J ' +1 ; the lowest 
mixed terms arc T n+1 t* and T(r*) j+1 . 

The above argument lays the foundation why the 
RUDD outperforms other sequences. To illustrate the ar- 
gument we plot the dependence of A p p on r* for various 
sequences of finite-width pulses in Fig. [5] for the central 
spin model at M = 8. Results for the spin chain model 
(not shown) look very much the same except for a rescal- 
ing of T. In panel (a) all sequences behave similarly; the 
dependence on r* is linear, and the RUDD behaves worst. 
This fact is attributed to the larger average length of the 
pulses. Note that in the regime depicted the distance 
A p F is still fully dominated by the pulse errors. 

In panel (b) we can nicely see the crossover from the 
regime where the pulse error dominates (straight lines 
corresponding to (t*) 2 ) to the saturation levels corre- 
sponding to the errors of the ideal CPMG and CDD se- 
quence. The errors of the ideal UDD sequence is much 
lower so that its saturation level cannot be seen. Still the 
TiRUDD behaves worse than the TiUDD. 
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FIG. 5: (Color online) Distance A p f vs. the shortest pulse 
duration r* for various sequences at a = 10 with T = 0.09/A; 
the panels correspond to pulses of finite width of different 
order as in Fig. [3] 



In panel (c) we again see the crossover from pulse errors 
to sequence errors on t* — > 0. Interestingly, the RUDD 
behaves better than the UDD in that the pulse errors 
decrease expectedly faster A p f,rudd °c ( t *) 3 compared 
to App, udd fx ( T *) 2 - We stress that the latter scaling is 
no contradiction to the pulse being second order because 
an error T(r*) 2 is not excluded. Fig.EJc) establishes that 
such mixed terms indeed deteriorate the performance of 
unadaptcd sequences of finite- width pulses. This clarifies 
the behavior of RUDD relative to other sequences. 

For practical implementation, it is important to point 
out that the behavior of T2RUDD for small r* is indepen- 
dent of whether or not we include the very short bound- 
ary 27T pulses, cf. solid line and circles in Fig. EJc). This 
is due to the shortness of these effective identity pulses. 

Last but not least, we find another regime of low val- 
ues of A P F- This is the regime where the pulse lengths 
reach their maximum value because the pulses touch one 
another. They are back to back. Quite unexpectedly, the 
full RUDD including the boundary 27r pulses again per- 
mits to obtain an extremely good suppression of decoher- 
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cncc. This regime is very interesting because it requires 
only very low pulse amplitudes and a small total energy 
for the coherent control, cf. Eq. (fTU|) . due to the pulses of 
maximum length. Further studies of this relevant regime 
are left to future research. 

Conclusions. The analysis of sequences of finite- 
width pulses allows us to draw the following conclusions. 
They are derived from the data for the models studied, 
but we expect them to hold more generally. 

First, the use of higher order pulses generically implies 
a significant improvement. Such pulses are designed such 
that they suppress the coupling to the bath to a high 



order during their action |25(. Second, non-equidistant 



sequences such as UDD outperform or, in the worst case, 
perform the same as equidistant (CPMG) or concate- 
nated (CDD) sequences. 

Third, in the regime, where the pulse errors dominate 
the suppression of decoherence is further enhanced by 
varying the pulse durations alon g th e sequence (RUDD) 
as suggested on analytic grounds [271 ] . This enhancement 
takes only place for pulses of sufficient high order. We 
found that it is present for second order pulses. This 
establishes RUDD as a promising concept and represents 
our central result. 

Fourth, an additional interesting asset of the RUDD 
is that the total energy required for the coherent control 
by pulses increases only logarithmically with the number 
of pulses - in contrast to all other sequences of unvaried 
pulses. Hence in particular long coherence times can be 
realized at low energy price. 

Fifth, surprisingly, we found an additional regime 
where the RUDD suppresses decoherence efficiently. This 
is the regime where the pulses are (almost) back-to-back 
approaching continuous modulation Q. Because in this 
regime the pulses reach their maximum length the re- 
quired control energy is a minimum. Further research is 
required to study this promising regime in detail. 
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